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Abstract 

The analytical bounce solution is derived in terms of the polygamma function 
in the Caldeira-Leggett's dissipative quantum tunneling model. The classical 
action for the bounce solution lies between the upper and lower bounds in 
the full range of a, where a is a dissipation coefficient. The bounce peak 
point increases from 1 to 4/3 with increase of a. In spite of various nice 
features we have shown that the solution we have derived is not exact one by 
making use of the zero mode argument in the linearized fluctuation equation. 
However, our solution can be a starting point for approximate computation 
of the prefactor in this model. 
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How to describe a dissipation at the level of quantum mechanics is a long-standing 
puzzle in physics. Upon our knowledge Feynman and Vernon(FV) [1] firstly described it as 
an interaction between a system of interest and its enviroments. Especially, they developed a 
formalism to investigate the quantum dissipation systematically by introducing an influence 
functional. The influence functional is an extremely important quantity in the sense that it 
contains all quantum effects of the enviroments and, thus makes it possible to describe the 
quantum dissipation in terms of only system's coordinates. 

Based on FV formalism Caldeira and Leggett(CL) considered in Ref. [2] a quantum tun- 
neling model interacting with the harmonic oscillator enviroments. Introducing a potential 
which does not have a true vacuum, CL examined the effect of dissipation on quantum tun- 
neling within a semi-classical approximation which was developed several decades ago [3-5] . 
The full application of the semi-classical application, however, is very difficult in this setup 
because the exact bounce solution is still unknown. In fact, it seems to be extremely hard 
(or may be impossible) to derive a bounce solution in an analytic form due to the non-local 
term the CL model involves. Without the analytic bounce it is impossible to exploit the 
power of the semi-classical approximation maximally. 

In Ref. [2] CL obtained the upper and lower bounds of the classical action in terms 
of the dissipation coefficient a without the explicit bounce solution. These two bounds 
are monotonically increasing function with respect to a, which indicates that the presence 
of the dissipation causes the decrease of the tunneling probability within an exponential 
approximation. The similar physical setup with a double-well potential was examined by 
a canonical method in Ref. [6,7]. The authors in Ref. [6,7], however, claimed that the 
dissipation may enhance the tunneling probability. In order to reconcile these two apparently 
discrepant results we think the prefactor should be examined in the semi-classical side. 
However, it is very difficult to compute the prefactor without the analytic bounce solution. 
Thus we may need a bounce solution in the analytic form for the computation of the prefactor 
although it is not exact. It is main purpose of this letter to derive an analytic bounce which 
interpolates between the exact no-damping solution and the strong-damping solution. 
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We start with a dimensionless action [2] 1 
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which yields an equation of motion as 
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We know the exact no-damping (a — > 0) solution 
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and strong-damping (ct — > oo) solution 
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Thus, the real bounce solution should interpolates between (3) and (4) with increasing the 
dissipation coefficient a. 

The bounce solution we obtained in this letter is 
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where ip' is an usual polygamma function. The A a , B a , and C a are a-dependent but u 
independent constants which obey the following equations: 
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(3 - 8a 2 )/iV(X) + 6(1 - 4a 2 )/iV(/u) + 16a 2 /i + 4a 2 = 

where ■0" is polygamma function and /i = C a /B a . From the first and second equations of 
Eq.(6) one can show easily C a = 2a. One can show also numerically that A a and B a are 
monotonically increasing functions with respect to a. 



In this letter we will follow the same conventions with Ref. [2] 
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Before we explain how Eq.(5) is derived, we would like to show its nice features. Fig. 1 
shows the a-dependence of its classical solution which is represented by the red line. The 
green and blue lines represent the classical actions for the strong-damping and no-damping 
solutions respectively. Two black lines are upper and lower bounds of the classical action 
which were derived explicitly in Ref. [2]. Fig. 1 indicates that the classical action for our 
bounce solution (5) lies between the upper and lower bounds in the full range of a. Fig. 1 
also implies that our solution interpolates between the no-damping and the strong-damping 
solutions. The small difference between the red and green lines at the large a region indicates 
that our solution is an approximate analytical solution. This fact will be proven later by 
making use of the zero mode argument in the fluctuation equation level. Fig. 2 shows the 
a-dependence of -2(0), which means the peak point of the bounce. As CL predicted in Ref. 
[2], the peak point increases from 1 to 4/3 with increase of a. 

Now, let me explain how the bounce solution (5) is derived. Taking a Fourier transform 
from z(u) to z{uj), one can change the equation of motion (2) in terms of z{uS) as following 



The explicit form of Zq(u) and z^uj) which are Fourier transform of z (u) and z^iu) become 



As expected z (uj) and z^uo) are solutions of Eq.(7) without second term and first term in 
the left-handed side of Eq.(7) respectively. These are properties of no-damping and strong- 
damping in cu-space. It is worthwhile noting the ^-dependence of zq{uj) and z^uj). Although 
uj is a dimensionless quantity, it should have a dimension if we go back to the CL's original 
model which has a dimension. Thus, from the dimensional consideration in Eq.(7) we can 
conjecture that z{uS) should have a same dimension with uj in the original theory except the 
strong-damping solution, which should be dimensionless. From this point of view we can 
understand the ^-dependence of zq(u) and z^uj). Thus one can take an ansatz 
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Then the a-dependent constants A a , B a , and C a should satisfy A = 2y/2~jr, B = n, C = 0, 
Aoc/Bcv = 4a/27Tq;/3 and = 2a in order for z(uj) to interpolate between the no-damping 
and the strong-damping solutions. 

Now, we will insert the ansatz (9) into Eq.(7) to extract an information on A a , B a , and 
C a . The most difficult term we need to compute is the following convolution term 
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Note that Ii ja (w) and l2, a (w) are even function with respect to uj. Thus we can assume 
uj > without loss of generality. With this assumption Ii tCC (u) and l2, a (w) are expressed as 
following 
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and z = i? Q cj and /i = C a /B a . When we compute the first term of Ii, a {uj) we used the 
property of the Lerch function in Ref . [8] . 

Now, the remaining problem for the computation of I a (u) is to compute K which have 
an infrared-like infinity as a field theory terminology. In order to take into account the 
infinity carefully we take a change of variable x = e z , which makes K to be 
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where xo = e z °. In Eq.(14) we introduced an infinitesimal parameter e explicitly for the 
regularization of the infrared-like infinity. Performing the integration in Eq.(14) one can 
express K as a difference of two hypergeometric functions. Making use of the relation 
between the hypergeometric and digamma function [9] the final expression of K becomes 
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where (//)„ = + 1) •••(// + n — 1) and ^ is a digamma function. Note that K has a 
logarithmic divergence as expected. Using Eq.(15) it is straightforward to compute l2, a (w) 
which reduces to 
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Note that the infinity term in Eq.(15) disappears in Eq.(16) because of the exact cancellation. 
This exact cancellation also takes place in I 1>a {u). After tedious calculation the final form 
of I a (u>) reduces to 
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where 7 is an Euler's constant. In order for z{uj) to be an exact solution the right-handed 
side of Eq.(17) should be equal to 2v / 2~7r(u; 2 + 2au + l)/3. To extract an information on A a , 
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B a , and C a we assume this equality. Repeating to take a u — > limit and subsequently to 
differentiate the right-handed side of Eq.(17) three times, we can derive Eq.(6). Taking an 
inverse- Fourier transform to z(cu), we can derive Eq.(5). 

Although our bounce solution has many nice features as discussed before, it is not an 
exact solution unfortunately except a = 0. For a = case we can show analytically 
Aq = 2y / 27r, B = 7r and C = by solving Eq.(6) in the /i — > limit. This means 
our solution exactly coincides with no-damping solution at a — > limit. However, for the 
nonzero a we can show that our bounce solution is an approximate one by using the zero 
mode of the linearized fluctuation equation as following. 

Note that the CL model has a time-translational symmetry in spite of the presence of the 

non-local term. Although one can show this fact simply from the equation of motion (2), we 

would like to show it at the fluctuation level for a later use. Inserting z(u) = z d (u) + r)(u) 

into Eq.(2) one can construct easily the linearized fluctuation equation 

2a f°° , .riiu) - n(u') , . , 
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where A is an eigenvalue of the fluctuation equation and the constants ci, c 2 , and C3 are 
introduced for convenience. If c\ — c 2 — 1 and C3 = 0, Eq.(18) is a fluctuation equation 
around the no-damping solution. In this case it is trivial to show that dz (u)/du is a zero 
mode. If Ci = and c 2 = C3 = 1, Eq.(18) corresponds to a fluctuation around the strong- 
damping solution. In this case also one can show directly that dz^uj/du is a zero mode. 
This means that the CL model has a time-translational symmetry in spite of the presence 
of the dissipation. 

Now, let us consider the full fluctuation equation, i.e. c± — c 2 = c 3 = 1. If our bounce 
solution (5) is an exact one, by same reason 
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should be a zero mode. Inserting Eq.(19) into the non-local term in Eq.(18) one can show 
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where ((p, q) is a Riemann Zeta function defined as ((p, q) = J2T=o V(9 + k) p - 

Using Eq.(20) the left-handed side of the fluctuation can be plotted numerically. When 
a is small, the left-handed side of the fluctuation equation is plotted in Fig. 3 in terms of 
u for various a, which indicates that i] is not exact zero mode although approximately it is. 
Fig. 3 also shows how rj can be a zero mode in the a — > limit. Fig. 4 is a plot of the 
left-handed side of Eq.(18) when a is large. This figure also shows how our solution goes to 
the strong-damping solution in the a — > oo limit. But we should comment that our bounce 
solution does not seem to exactly coincide with the strong-damping solution in the a — > oo 
limit. That is why there is a small difference in Fig.l between the classical action for our 
solution and its lower bound at the large a regime. 

In this letter we have derived the analytic bounce solution in the CL model. The classical 
action for our solution lies between the upper and lower bounds in the full range of a. 
Although it has many nice features, we have shown that it is not an exact solution except 
a = case. However, using this approximate solution, one may be able to compute the 
prefactor approximately. We guess this prefactor may be important factor to reconcile the 
discrepancy between the semi-classical method and the canonocal method. We hope to 
visit this issue in the near future. Our analytic bounce solution might be extended to the 
non-Ohmic dissipation case. The explicit result will be discussed elsewhere in detail. 
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FIGURES 

FIG. 1. Plot of a-dependence of classical actions for our bounce solution (5) (red line), 
no-damping solution (3) (blue line) and strong-damping solution (4) (green line). Two black lines 
are upper and lower bounds of the classical action which were derived by CL in Ref.[2] without an 
explicit bounce solution. This figure indicates that the classical action for our solution lies between 
the upper and lower limits in the full range of a. 

FIG. 2. Plot of a-dependence of the bounce peak point. As CL have argued in Ref.[2], the 
bounce peak point increases from 1 to 4/3 with increase of a. 

FIG. 3. The plot of the left-handed side of Eq.(18) for various small a. The difference from 
zero indicates that our solution is not an exact one. The increase of the peak point and the 
wide-spreading shape of peak with increase of a denotes that our bounce solution (5) goes away 
from the exact one when a becomes larger in the small a regime. 

FIG. 4. Plot of the left-handed side of Eq.(18) for various large a. The decrease of the peak 
point with increase of a means that our solution approaches to the exact one when a becomes 
larger in the large a regime. 
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